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ABSTRACT 

We describe the High Precision Polarimetric Instrument (HIPPI), a polarimeter built 
at UNSW Australia and used on the Anglo-Australian Telescope (AAT). HIPPI is an 
aperture polarimeter using a ferro-electric liquid crystal modulator. HIPPI measures 
the linear polarization of starlight with a sensitivity in fractional polarization of ~4 
x 1CU 6 on low polarization objects and a precision of better than 0.01% on highly 
polarized stars. The detectors have a high dynamic range allowing observations of the 
brightest stars in the sky as well as much fainter objects. The telescope polarization 
of the AAT is found to be 48 ± 5 x 10” 6 in the g' band. 

Key words: polarization - instrumentation; polarimeters - techniques: polarimetric. 


1 INTRODUCTION 


Stellar polarization can be measured to very high sensitivity 
with ground-based telescopes. Unlike photometry where at¬ 
mospheric effects limit the achievable precision, polarimetry 
is a differential measurem ent and so there is no fundamen¬ 
tal limit on the sensitivity. iKemp et al.l d 1987tl measured the 
polarization of the Sun to levels of parts in ten million, and 
astronomical polarimeters have been built that measure stel ¬ 
lar polarization at the parts per mi llion level (I Ho ugh et al.l 
1200 j ; IWiktorowicz fc Matthewsl 120081 1. Interest in such in¬ 
struments has been driven, in particular, by the possibility 
of detecting polarized scattered light from hot Jupiter type 
exoplanets, which is predicted to be at levels of ~ 10~ 5 
or less in the combined light of the star and the planet 
llSeager. Whitney fc Sasselovll2000l : iLucas et alJl2009h . Such 


instruments also have other ap plications such as the study of 
the local interstel lar medium (|Bailev, Lucas fc Houghll2010l : 
IPrisclmjKd. 2012) a nd th e scattered light from debris disks 
I Wiktorowicz et al. HmS). 

The standard method used for most high sensitivity po¬ 
larization studies has been the use of the photoelastic mod¬ 
ulator (PEM) technology pioneered by James K emp llKemd 
Il969l : IKemp fc Barbounll98ll : IKemp et al .111 9871 ). These de¬ 
vices modulate polarization at frequencies of 20 - 100 kHz 
by means of the stress birefringence in an optical material 
made to vibrate at its natural frequency using piezoelec- 
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trie transducers. While PEMs have proved very successful in 
this role they nevertheless have a number of disadvantages. 
While a high modulation frequency is desirable to minimize 
the effects of variations due to seeing and tracking errors, 
the frequency of PEMs is higher than is really needed for 
this purpose. The high frequency can present problems in 
providing a suitable detector system. Many detector types 
can either not operate at the required speed, or can only do 
so with some compromise in their noise performance. PEMs 
are inherently sine-wave modulators and thus the efficiency 
of a PEM polarimeter is less than that of an ideal square 
wave modulator by a factor of at least y/2 (actually slightly 
more than this, see lHough et al.l l 200 6‘). PEMs are also quite 
bulky and thus difficult to use where space is limited. 

In this paper we describe a high-sensitivity polarime¬ 
ter using an alternate type of polarization modulator, a 
ferro-electric liquid crystal modulator (FLC). FLCs are elec¬ 
trically switchable wave plates consisting of a thin layer 
of liquid crystal material sandwiched between two glass 
plates. They have a fixed retardation and the orientation 
of the optical axis can be controlled by an applied drive 
voltage. FLCs have been used , in particular, for solar po¬ 
larimetry (e.g. iGandorferlIl999l : iMartmez Pillet et al.l 1 19991 : 
Hanaokj|2004|) and als o for the Extreme Polarimeter (ExPo 
Ro denhu is et al. [mil), a high contrast imaging polarimeter 
and the SPH ERE/ ZIMPOL instrument for the ESO VLT 
iBazzon et al .1 [20 1 j iRoelfsema et alJl2014h . 

Our polarimeter HIPPI (High Precision Polarimetric 
Instrument) has been used on the 3.9m Anglo-Australian 
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wheel also includes a clear position and a blank setting that 
can be used for taking dark measurements. 

The polarization analyser is a calcite Wollaston prism 
that provides a 20 degree beam separation. This is placed 
between two lenses to collimate the light through the prism. 
A Fabry lens in each beam images the telescope pupil onto 
the two detectors. The whole optical system from the colli¬ 
mating lens to the detectors is rotatable about the optical 
axis using a Thorlabs NR360S NanoRotator stage. Rotating 
this system through 90 degrees relative to the modulator 
has the effect of reversing the sign of the modulation seen 
by the detectors and provides a “second stage chopping” 
which helps to improve accuracy by eli minating some sys¬ 
tematic effects dKemp fcj3arbour 198llh A similar system 
was used in PlanetPol ( Hough et al.ll2006h . All the optics 
are anti-reflection coated for the wavelength range 350 
700 nm. 

The instrument components are mounted on a standard 
300 mm square aluminium optical breadboard that is at¬ 
tached by 90 degree brackets to a mounting plate that bolts 
to the back of the telescope. Many of the structural compo¬ 
nents, optical mounts and electronics enclosures have been 
constructed by 3D printing in ABS plastic. The instrument 
is therefore compact and lightweight (10 kg). 


Figure 1 . Schematic diagram of HIPPI optical system (not to 
scale) 


Telescope (AAT) at Siding Spring Observatory for two 
observing runs during 2014. The design of HIPPI has 
been based o n that of the previous PlanetPol instrument 
flHough et alJ [200~6 L However, HIPPI differs from Planet¬ 
Pol in being optimised for observation at blue wavel e ngths . 
Thi s is based on result s such as those of I Pont et al.l (l2013l l 
and I Evans et alJ d2013l l that indicate the presence of strong 
Rayleigh scattering at blue wavelengths from the exoplanet 
HD 189733b making these wavelengths the most suitable for 
detecting exoplanet polarization. 

2 INSTRUMENT DESCRIPTION 
2.1 Overview 

A schematic diagram of the HIPPI optical system is given 
in figure [T] The FLC modulator is the first element in the 
optical system. This is an important design feature since 
any optics placed ahead of the modulator could potentially 
induce spurious polarization effects, for example, polariza¬ 
tion due to inclined mirrors or residual stress birefringence 
in refracting elements. Following the FLC is an aperture of 
1 mm diameter corresponding to 6.7 arc seconds at the AAT 
f/8 focus. This is followed by a six position filter wheel. 

The filters used with HIPPI have been Sloan Digital Sky 
Survey (SPSS. iFukueita et alJll996h g' and r' filters (from 
Omega Optical, giving wavelength ranges of ~400 - 550 nm 
and ~550 - 700 nm respectively). There is also a short pass 
filter which passes wavelengths shorter than 500 nm (re¬ 
ferred to as 500SP). The instrument has little throughput 
below about 350 nm due to absoprtion in the calcite prism 
so the range of this filter is from ~350 - 500 nm. The filter 


2.2 Ferro-electric liquid crystal modulators 

Two different FLC modulators have been used with HIPPI. 
The first is a LV1300-AR-OEM device from Micron Technol¬ 
ogy it is designed for the 400-700 nm range and is 12.7 mm 
in diameter housed in a 25 mm diameter cell. The second is 
an MS Series polarization rotator from Boulder Nonlinear 
Systems (BNS) designed for the wavelength range 425 - 675 
nm and is 22 mm in diameter with a 15 mm useful aper¬ 
ture. Both devices are designed to be half-wave retarders at 
a wavelength near 500 nm, and depart from half-wave away 
from this wavelength as discussed further in section [3761 

The two modulators are very similar in their operation 
and provide good polarization modulation with a ±5 V drive 
waveform. However, we have found the BNS modulator to 
produce much lower levels of instrumental polarization, and 
it is therefore currently the preferred option. 

Electrically the modulators are equivalent to capacitors 
of ~200 nF and therefore require a drive circuit capable of 
driving at high speed into a capacitive load. The devices can 
also be damaged by sustained DC voltages. We have desiged 
and built a drive circuit consisting of a two-pole Butterworth 
high pass filter followed by an amplifier using a NPN/PNP 
transistor pair output stage. The filter ensures no DC or low 
frequency components reach the device. The drive amplifier 
has the high slew rate, and high drive current needed to 
drive a square wave into the capacitive load. 

The drive waveforms are generated in software. A sim¬ 
ple square wave between +5 V and —5 V has been used 
for all the observations described in this paper. Our system 
allows selection of modulation frequencies between 200 Hz 
and 2 kHz. We have found 500 Hz to be a good choice for 
actual observing, providing a close to square wave modula- 
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tion, while being fast enough to be insensitive to intensity 
fluctuations due to seeing or tracking errors. 

FLCs are temperature sensitive devices. The switching 
is faster at higher temperatures and the switching angle is 
also temperature dependent. To ensure consistent and stable 
operation we mount the FLC in a temperature controlled 
lens tube and operate it at a constant temperaure of 25° C, 
maintained to about ±0.1° C. 


2.3 Detectors 

The detectors used in HIPPI are compact photomultiplier 
tube (PMT) modules. The modules contain a metal pack¬ 
aged photomultiplier tube combined with an integrated FIT 
supply. PMTs have substantial advantages of large detector 
area and low dark noise compared with possible solid-state 
alternatives such as avalanche photodiodes (as used in Plan- 
etPol) or so-called “silicon photomultipliers”. 

FIIPPI uses Flamamatsu F110720-210 PMT modules 
which have ultra-bialkali photocathodes dNakamura et al] 
l2010h providing a quantum efficiency (QE) of 43% at 400 
nm. They are compact modules operating off a single 5V 
supply. For the high photon rates required with HIPPI it is 
not possible to use the PMT in a photon counting mode. 
Instead we use a transimpedance amplifier to amplify the 
photo-current. The amplifiers designed and built for HIPPI 
use an ultra low noise Texas Instruments OPA 129 opera¬ 
tional amplifier with an extremely low input current noise of 
0.1 fA Hz -1 / 2 . Remotely switchable transimpedance gains 
of 10 5 , 10® and 10 7 V/A can be selected. The PMT module 
itself provides selectable HT voltages from 500 V to 1100 V 
corresponding to a variation of the PMT gain from 5 x 10 3 
to 3 x 10® electrons per photon. The ability to remotely 
vary both the photomultiplier gain and amplifier gain over 
a wide range provides a very high dynamic range, allowing 
HIPPI to observe objects from the brightest stars in the sky 
to quite faint objects while still providing close to photon 
noise limited performance. 

The PMT amplifiers have been built in surface-mount 
construction on a compact printed circuit board 25 x 50 
mm that fits on the back of the PMT module as shown in 
figure [5] 


2.4 Instrument Control and Data Acquisition 

HIPPI is controlled by software running on a rack mount 
computer (Intel quad core i7, 8 GB RAM, 2 1 TB disks). 
The computer runs the Windows 7 Professional operating 
system and the software has been developed using the Na¬ 
tional Instruments Lab VIEW graphical programming envi¬ 
ronment. 

The interface to drive the FLC and read data from the 
detectors makes use of two National Instruments data ac¬ 
quisition modules (PCIe 6341) each of which provides 16-bit 
analog input and output channels. The drive waveform for 
the modulator is generated in software and output from one 
of the modules. A trigger digital input signal is generated 
from the rising edge of the square wave and fed to both mod¬ 
ules. The detector signals are read by analog input channels 
on each module. A schematic diagram of the system is shown 
in figure [3] 



Figure 2. HIPPI Detector Modules 


The HIPPI software system also provides control of the 
FLC temperature, the filter wheel selection, the rotation of 
the Wollaston prism and detector section of the optics, and 
the gain and HT voltage settings for the detectors. 

In operation the analog input channels are sampled at 
10 microsecond intervals and read for an integration time of 
typically 1 second resulting in 100000 data points for each 
channel. The timing is controlled to start on the trigger in¬ 
put so that the sampling always has the same phase relation¬ 
ship with the modulator waveform. The data are then folded 
in software over the modulation cycle. For our standard 500 
Hz modulation this results in an array of 200 points from 
each channel. 

Figure [4] shows what the resulting waveforms look like 
when a polarized star is being observed. The modulation is 
of opposite sign in the two channels as they correspond to 
the two orthogonal polarization states from the Wollaston 
prism. Zero on the diagram corresponds to the rising edge 
of the square wave drive signal to the FLC. The delay of 
about 100 /rs in the observed signal is a combination of the 
finite switching time of the FLC, and the time constant (22 
/rs) of the detector amplifiers. 

These modulation waveforms are displayed to the user 
and output to the data files every integration. The ampli¬ 
tude of the modulation is a measure of the polarization of 
the source. A quick-look data reduction capability is built 
into the observing software and derives the fractional polar¬ 
ization (p) as follows. 


(A + Y) 


( 1 ) 


where X is the signal from the first flat part of the wave¬ 
form (points from about 20 100 in figure [4]), and Y is the 

signal from the second flat section (about 120 - 200), with 
points around the two transitions being ignored. This simple 
quick-look reduction allows immediate evaluation of incom¬ 
ing data, but excludes some corrections such as subtracting 
the small bias and dark signals. The offline data reduction 
described in section [3] uses the full waveform and inlcudes 
all these corrections. 
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Figure 3. Schematic diagram of HIPPI data acquistion system 
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Figure 4. Example of the modulated signal seen by HIPPI when 
observing a polarized star (in this case HD 147084 with about 
3.2% polarization). This is the observed signal integrated over 
one second and folded over the 500 Hz modulation cycle giving 
200 points over the modulation cycle. The rms noise on each point 
is 2.1 mV giving a S/N of 850. 


0 50 100 150 

Time (s) 

Figure 5. Second stage chopping procedure showing the sign of 
the measured modulation reversing for each 90 degree rotation 
of the Wollaston prism and detectors. The top panel is for the 
highly polarized star HD 147084 (about 3.2% polarization in this 
case), and the lower panel is for a lower polarizaton object (about 
150 ppm, or 0.015 %). 
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2.5 Observing Procedure 

The normal observing procedure with HIPPI is to make a 
sequence with the Wollaston and detectors rotated to two 
positions 90 degrees apart (referred to as A and B), in the 
cycle A, B, B, A with this sequence then being repeated as 
many times as required. This rotation causes the sign of the 
modulation in each channel to reverse. It results in the p 
value obtained for each integration (with Equation [T] above) 
following the sequence as shown in figure [5] 

Note that while the polarization changes are almost 
symmetric about zero for the highly polarized star in the 
top panel, a significant offset from zero is apparent in the 
bottom panel. 

An observation of this type measures only one Stokes 
parameter of linear polarization. To get the orthogonal 
Stokes parameter the entire instrument is rotated through 
45 degrees and the sequence repeated. The rotation is per¬ 
formed using the AAT’s Cassegrain instrument rotator. In 
practice we also repeat the observations at telescope position 
angles of 90 and 135 degrees to allow removal of instrumen¬ 
tal polarization. 


3 DATA REDUCTION 

The HIPPI data processing and calibration procedure has 
at its heart a Mueller matrix model of the instrument. The 
parameters of the Mueller matrix are varied according to 
the state of the system. The system parameters vary over a 
single integration as a result of the varying voltage applied 
to the FLC. 

The Mueller matrix M of an optical system (such as 
HIPPI) relates the output Stokes vector s ou t to the input 
Stokes vector s; n through: 


Sout = Ms in (2) 

The Mueller matrix for HIPPI is not constant since it 
changes through the cycle of the modulator. We can describe 
a full integration using a system matrix W. The system ma¬ 
trix is an N by 4 matrix, where each row is a state of the 
system, corresponding to a single data point in the mod¬ 
ulation curve. Multiplying the input Stokes vector by the 
system matrix gives the vector x of N observed intensities 
seen at the detector during the modulation cycle (e.g. the 
N points plotted in figure [4] where N = 200). 


X = Ws in (3) 

It can be seen that the N rows of the system matrix 
W are simply the top rows of the Mueller matrices M for 
HIPPI for each of the N states of the instrument through a 
modulation cycle. Only the top row is needed because this is 
the row of the Mueller matrix that determines the intensity 
of the output light, and only intensity is directly measured 
at the detector. 

The optical components HIPPI is made up of are the 
FLC and the Wollaston prism, the FLC is a retarder, which 


has the Mueller matrix: 


0 0 
5(7(1 — cos 5) — 5 sin 8 

S 2 + C 2 cos 5 C sin 5 
—(7 sin 5 cos 5 

(4) 

where C = cos 20, 5 = sin 20, S is the retardance, and 0 the 
angle to the fast axis of the retarder. For a half wave plate 
the retardance is n radians. FLCs are someti mes mo delled 
as h aving a linear depolarization component iGendre et ali 
l2010ll . with a Mueller matrix of: 


1 0 
0 C 2 + S 2 cos <5 
0 5(7(1 — cos 5) 

0 5 sin 5 


Mloepol 


10 0 0 

0 1 — d 0 0 

0 0 1 -d 0 

0 0 0 1 — d 


(5) 


where d is the degree of depolarization. 

The Wollaston prism behaves as two perpendicular po¬ 
larizers. A polarizer has a Mueller matrix of: 


1 


Mpoi = 1/2 


eC 

eS 


0 


eC 

eC 2 

eCS 

0 


eS 

eSC 

eS 2 

0 


0 

0 

0 

0 


( 6 ) 


where (7 = cos2i9, 5 = sin 2d, e is the efficiency of the 
polarizer, and d is the angle of the polarizer axis to that 
defined for the incoming beam. Wollaston prisms are very 
efficient, and we assume e = 1. 

The combined Mueller matrix for HIPPI can now be 
written in terms of these matrices as: 


M = MpolMRetMQepol (7) 

and can be used to derive the system matrix for HIPPI as 
already described. The N rows of the system matrix are de¬ 
rived from the Mueller matrices for the state of the system 
at each modulation point where each row differs due to dif¬ 
ferent values of 0 (in Mpet) and d (in Mpepoi). 

The other main parameter in the Mueller matrix is the 
retardance 5. This does not vary around the modulation 
cycle, but it is a strong function of wavelength. The re¬ 
tardance is close to half-wave only at a single wavelength. 
The effect of a retardance that is not exactly half-wave is to 
reduce the amplitude of the modulation curve by a factor 
of (1 — cos<5)/2. Incorporating the full wavelength depen¬ 
dence of retardance in our reduction model would be im¬ 
practicable as a single observation can cover a wide range of 
wavelengths. Instead we perform the reduction assuming a 
half-wave retardance, and then account for the wavelength 
dependence effects by scaling the polarization with an effi¬ 
ciency correction that can be determined empirically from 
standard star observations (section 14.311 or can be derived 
from a model ('section |3.6I) . 

On ce t h e system matri x is knwn equation [3] can be in¬ 
verted (ISabatke et al.ll2000l ) to give: 


s in = W+x (8) 

where W + is the pseudo-inverse of W, which 
we c alculate numeri cally using the method of 
iRust, Burrus fe Schneeburgeii (Il966li . This gives the 
source Stokes parameters Si n in terms of the observed 
modulation data x. 
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Figure 6. The structure of a measurement 


Ill practice, because HIPPI isn’t designed to mea¬ 
sure circular polarization, we use only the first three 
columns of W + , to obtain a least squares estimate of 
th e components o f Si n , ( I, Q and U), in the manner of 
lOendre. Foulonneau fc Bigud 1 20101 1 and then divide Q and 
U by I for the normalised Stokes parameters. 


3.1 Calibration 

To apply equation[ 8 ]we need to know how the waveplate an¬ 
gle (j) and the depolarization factor d vary thrugh the mod¬ 
ulation cycle of the FLC as these are the only varying com¬ 
ponents that enter into W. This is done using a laboratory 
calibration procedure where we input known polarization 
states into the instrument using a lamp and a rotatable po¬ 
larizer. Measurements are made of the full modulation curve 
with the polarizer stepped through a full rotation in 10 or 20 
degree steps. By fitting a Malus law (jHec htlteoOlh — mod¬ 
ified to include an intensity offset and a phase shift — to 
the intensity as a function of polarizer angle we determine 
</> and d for each point in the modulation cycle. 

To reduce the effect of noise in the calibration measure¬ 
ments the array of cj> values is smoothed by replacing the last 
50 modulation points in the plateau regions (i.e. the last 50 
before the FLC is switched) by their average value. Only the 
last half of the plateau region is used because we have found 
that the Micron Technology FLC does not have a stable <j> 
until this point - it is still approaching maximum/minimum. 

The values of cj> and d are determined separately for 
each of the two PMT channels. This is to allow for timing 
differences between the two channels. The depolarization (d) 
is included primarily to account for the switching part of the 
modulation cycle, where the switching can occur faster than 
the detector response resulting in a reduction in apparent 
polarization. The results are normalized such that it does 
not result in any scaling of the measured polarization, and 
therefore does not attempt to duplicate the wavelength de¬ 
pendent efficiency correction. 


3.2 Measurements 

For a given filter and target, a measurement involves taking 
data at 4 position angles of the Cassegrain rotator corre¬ 
sponding to 0, 45, 90 and 135 degrees. In general, a sky 


measurement is also made at each telescope position an¬ 
gle. For very bright targets or high polarization standards 
or bright targets in moonless conditions a dark measure¬ 
ment may be substituted for the sky measurements. Each 
sequence is made up of a number of identical repeats each 
of which consist of a subsequence of a number of rotations; 
for each rotator position there are a number of integrations, 
each of which consists of the individual modulation points 
that make up the modulation curve. The rotation subse¬ 
quence is typically an ABBA sequence with A and B cor¬ 
responding to orthogonal rotator positions as described in 
section [53] Figure [S] summarises the components of a mea¬ 
surement. 

As an example, a typical bright star measurement with 
HIPPI consists of sequences at 4 position angles each con¬ 
sisting of 2 repeats (n rep = 2), in an ABBA rotation subse¬ 
quence (UrotA = UrotB = 2 ), each consisting of 10 1 second 
integrations ( m n t = 10 ) of a modulation curve consisting of 
200 modulation points ( n p t = 200). Measurements requiring 
higher precision have higher numbers of integrations and re¬ 
peats. 


3.3 Statistical treatment of values and errors 

The Mueller matrix model of the modulation curve is used 
to calculate a Q/I and U/I for each integration. However, 
the co-ordinate frame of the instrument is chosen to match 
the calibrated rotator positions (see section IXTI which means 
that the majority of the modulation points in the curve mea¬ 
sure predominantly either Q or U; the other Stokes parame¬ 
ter determination is discarded because it contributes greater 
noise. The “off-axis” Stokes parameter is measured by sig¬ 
nificantly fewer modulation points, and is contaminated by 
the intrinsic FLC polarization. This gives a single Stokes 
parameter determination for each integration. 

The average of all the integrations for a given rota¬ 
tor position is then calculated, and the error calculated 
from the standard deviation of the m„t points. The aver¬ 
age polarization over each rotation, repeat and channel in 
the sequence is calculated, and the error on this calculated 
by averaging the statistical errors for each individual rota¬ 
tor position and divding by the square root of the num¬ 
ber of these that have been combined. This gives a nor¬ 
malised Stokes parameter measurement for each telescope 
position angle (So, S 45 , Sgo, S 135 ), each with an associated 
error (Eo, E 45 , Ego, E 135 ). 

The value of Q/I, in the chosen co-ordinate system of 
the instrument, is obtained from the average of the 0 and 
90 degree sequences and the value of U // from the average 
of the 45 and 135 degree sequences. 


Q/I = 


So + 5*90 


U/I = 


S 45 + Sl35 


With errors calculated as: 


Eq/i - 


E 0 + Ego 

2y/2 


(9) 

( 10 ) 


( 11 ) 
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, _ Ei 5 + E 135 

' U/I - 2V2 


( 12 ) 


3.4 Sky and dark correction 

Sky measurements are made to take account, predominantly, 
of polarization caused by reflected light from the Moon. The 
sky sequences are always made with the same HT voltage, 
gain, telescope position angle and filter as the science se¬ 
quences. Because the sky signal level is much lower than 
that from the star, much shorter total integration times can 
be used for the sky measurements without intoducing sig¬ 
nificant additional noise. The average modulation curves for 
the sky data are determined for each PMT channel and rota¬ 
tor position. These are then subtracted from the each of the 
corrsesponding modulation curves in the science sequence 
prior to the analysis described in sections 13.21 and 13.31 

In cases where the sky contribution is negligible, for ex¬ 
ample when observing a bright star in moonless conditions, 
a dark observation can be used as an alternative to a sky. 
This is made by using a blank in the filter wheel to exclude 
any light input to the instrument. 


3.5 Coordinate transformation and efficiency 
calibration 

The co-ordinate transformation from the frame of the instru¬ 
ment to that of the sky is carried out by reference to known 
polarization standards. Polarization angle is calculated by: 

0 = 1/2 arctan (13) 

'°G /-* 

The difference between the calculated angle with that 
known for the standards becomes the angular correction, 
and all other measurements are rotated by this value. 

The linear polarization, P, is given by: 

p = V(Q/i ) 2 + ( u/iy ( 14 ) 

This polarization then needs to be scaled by multiplying 
by an efficiency correction to account for the wavelength de¬ 
pendence of the retardance of the modulator. This efficiency 
correction can either be determined empirically using polar¬ 
ized standard stars as described in section POl or determined 
using the bandpass model described in section 13.61 


3.6 Filters and bandpass model 

The filters used in HIPPI are relatively broad (150 nm for 
the g' and r' filters) and hence the precise effective wave¬ 
length will vary with star colour and other factors, and this 
will, in turn, affect the modulation efficiency which is also a 
function of wavelength since the modulator is not an achro¬ 
matic device. To account for these effects we use a bandpass 
model similar to that described bv lHough et alj (l2006l l . 

As a sta rting point for su c h a model we typically 
use one of the lCastelli fe Kurucd d2004l f stellar atmosphere 
models. The spectral energy distribution (SED) given by 
the model can then be modifi ed for interstel lar extinction 
using the empirical model of ICardelli. Clayton fe Mathis! 
(j 1989T ). We then correct the SED for its passage through the 


Table 2. Properties of FLC Modulators 


Modulator 

Micron 

BNS 

Ao 

505 =b 5 nm 

498 + 5 rim 

Cd 

1.75 ± 0.05 x 10 7 nm 3 

1.70 ± 0.05 x 10 7 nm 3 


Earth’s atmosphere by applying an atmospheric transmis- 
sion correction calculated usin g the VSTAR modelling code 
({Bailev fc Kedziora-Chudczerl [2012 1 and including molecu¬ 
lar absorption and Rayleigh scattering. 

Finally this can be combined with a model of the instru¬ 
ment response which includes the transmission functions of 
the filters, and the cathode radiant sensitivity (in mA/W) 
of the PMT as taken from the Hamamatsu data sheet. The 
final result, which we call S( A) is the relative contribution 
to the output detector signal as a function of wavelength. 

The effective wavelength of the observation can be de¬ 
termined from 


A / XS{X)dX 
eff J S(X)dX 


(15) 


where the integral is taken over all wavelengths for 
which S(X) is non-zero. 

The polarization modulation efficiency correction will 
be given by 


fe(X)S(X)dX 
I S(X)dX 


(16) 


where e(A) is the wavelength dependence of the mod¬ 
ulation efficiency. For an FLC the efficiency as a func¬ 
tion of wavelength is primarily determined by the vari- 
ation of the retardance with wavelength. According to 
iGisler, Feller fe Gandorfei] (l2003l l the optical path difference 
(A) of an FLC can be represented by: 


a -t +cj 0-£) < i7 > 

where Ao is the half-wave wavelength, C is a parameter 
describing the dispersion in birefringence of the FLC ma¬ 
terial, and d is the thickness of the FLC layer. In practice 
Cd can be treated as a single parameter. We have deter¬ 
mined these parameters for our two modulators by fitting 
this formula to a set of laboratory measurements of a po¬ 
larized source made with HIPPI using a set of narrow band 
filters. The results are given in table [2] 

The properties of the two modulators are very similar, 
an d the Cd value we obtain is almost the same as that found 
bv lGisler et al.1 (2003) for a similar device. 

The modulation efficiency is then given by: 


e(A)~^=p(l —cos2 t4 ) (18) 

where e m ax is the peak efficiency measured at wave¬ 
length Ao- Ideally e m ax should equal 1, but we find a value 
of 0.98 is the maximum achievable using laboratory mea¬ 
surement with HIPPI. 
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Table 1. Effective wavelength and modulation efficiency for different spectral types according to bandpass model 


Spectral Type Effective Wavelength (nm) Modulation Efficiency (%) 



500SP 

g' 

r' 

500SP 

g' 

r' 

B0 V 

433.2 

472.6 

600.8 

74.3 

90.8 

84.4 

A0 V 

443.7 

475.6 

601.9 

81.6 

91.5 

84.2 

F0 V 

447.1 

479.5 

603.7 

83.0 

92.3 

83.8 

GO V 

451.2 

483.4 

605.5 

84.6 

93.0 

83.5 

K0 V 

456.4 

486.4 

606.6 

87.1 

93.7 

83.3 

M0 V 

459.2 

490.1 

610.5 

88.8 

93.8 

82.6 

M5 V 

458.1 

490.0 

611.4 

88.3 

93.7 

82.4 




--- - -, -— > 


- R Filter 

: *- 1 \ 


- - G Filter 

: ^ 


500SP Filter 

;• _ * ^ 


Clear 


o.o 



Wavelength (nm) 


Figure 7. Illustration of the various components that combine to give the effective filter bandpass (S(X)) shown in the lower panel of 
the plot 


Tabic [l] shows how the effective wavelength and po¬ 
larization efficiency vary with spectral type. The bandpass 
model is illustated in figure [7] 


4 PERFORMANCE AND RESULTS 

The performance of HIPPI has been evaluated based on our 
observing runs on the AAT over May 8 - 12 2014 and Aug 
28 - Sep 2 2014. 


4.1 Instrumental polarization 

We have found that the FLC modulators used in HIPPI 
introduce a large instrumental polarization effect which 
appears to be an intrinsic property of the modulator. 
Such polarization is thought to arise from multiple in¬ 
ternal reflections in the birefringent material between the 
plates resulting in wavelength dependent fringe patterns 
in transmittance that are different for l ight polarized par¬ 
allel or perpendicula r to t he fast axis llGisler et al.l 120031 : 
Ide Juan Ovelar et alj l2012l l. This effect was found to be 
present with both our modulators but is a factor of about 
3 higher for the Micron Technology modulator (used for the 


© 2015 RAS, MNRAS 000.HHT21 






















































A high-sensitivity polarimeter 9 


Table 3. Low polarization star measurements to determine tele¬ 
scope polarization 


Star 

Date 

Filter 

P (ppm) 

6 (deg) 

BS 5854 

Aug 28 

g’ 

35 ± 5 

113 ± 6 


Aug 29 

g' 

46 ± 4 

115 ± 17 


Aug 30 

g' 

50 ± 4 

114 ± 4 


Sep 2 

g' 

45 ± 5 

110 ± 6 


Average 

g' 

44 ± 2 

114 ± 3 

Beta Hyi 

Aug 28 

g' 

52 ± 3 

113 ± 4 


Aug 29 

g' 

62 ± 3 

106 ± 3 


Aug 30 

g' 

58 ± 3 

109 ± 4 


Aug 31 

g' 

52 ± 3 

106 ± 4 


Average 

g' 

56 ± 2 

109 ± 2 

Sirius 

Aug 31 

g' 

55 ± 1 

112 ± 1 


Sep 2 

g' 

48 ± 4 

108 ± 5 


Average 

g' 

51 ± 2 

110 ± 3 

Adopted TP 


g' 

48 ± 5 

111 ± 2 


May 8-12 observing run) than for the BNS modulator 
(used for the Aug 28 - Sep 2 run). Even in the BNS case 
the effect is at the ~1000 ppm level and highly variable with 
wavelength. Since it is due to the modulator that does not 
rotate, it cannot be removed by our second-stage chopping 
procedure described earlier. 

Fortunately because the instrument only measures one 
Stokes parameter at a time it is possible to arrange things 
such that the instrumental polarization is orthogonal to the 
Stokes parameter being measured. This is done by choos¬ 
ing optimal values for the pair of 90 degree separated angles 
that are used for rotating the Wollaston prism and detectors. 
When this is done the instrumental polarization is greatly 
reduced but it cannot be removed completely. Due to ef¬ 
fects such as its wavelength dependence in angle, residual 
instrumental effects at around the 50 ppm level remain. 

The residual instrumental polarization is removed by 
the procedure of observing at four position angles of the 
instrument rotator (0, 45, 90, 135) as described in sections 
I2.5l and l3.2l The 90 degree rotation between the 0, 90 and 45, 
135 pairs reverses the sign of the instrumental polarization 
relative to that from the star, enabling it to removed. 


4.2 Telescope polarization 

The telescope optics can be expected to introduce a small 
telescope polarization (TP) that must be corrected for. In 
the case of PlanetPol on the Wi lliam He rschel Tele scope TP 
was found to be around 15 ppm (iHough et al.ll2006l) . A much 
larger TP of ~ 250 ppm was found using POLISH on the 
5-m Hale Telescope (IWiktorowicz fc Matthewsli2008i ). 

As the AAT is an equatorially mounted telescope it is 
not possible to separate the t elescop e and star polarization 
in the way described bv lHough et al.l (l2006 f ) making use of 
the held rotation in an altazimuth mounted telescope. 

Instead, we have adopted as our preliminary estimate 
of the telescope polarization the average of our observations 
of three stars which we have reason to believe should have 
very low polarization. One of these is BS 5854 which is one 
of the low polarization standard stars observed with Plan¬ 


etPol. It is at a distance of 22.5 pc and had measured po¬ 
larizations with PlanetP ol of 5 ppm or lower (I Hough et al.l 
l200fil : iBailev et al.1 l2008l l . The other two are nearby bright 
stars Beta Hydri (HIP 2021, BS 98) which is a 2.8 magnitude 
star at a distance of 7.5 pc, and Sirius (a CMa) at a distance 
of 243 pc. Based on the polarization with distance results of 
IBailev. Lucas fc Houghl (|201dl l these should be expected to 
have very low polarizations. 

The results are given in table [3] for the telescope po¬ 
larization measurements in the SDSS g' filter. From smaller 
numbers of observations we have determined TP = 53 ± 2 
ppm, 9 = 110 ± 2 degrees in the 500SP filter, TP = 45 ± 
3, ppm 9 = 125 ± 5 degrees in the SDSS r' filter, and TP 
= 49 ± 2 ppm, 9 = 117 ± 2 degrees with no filter (the full 
wavelength range of the detector). These results suggest an 
increasing TP and a rotation in angle from red to blue. 


4.3 Polarized Standard Stars 

Observations were made of a number of polarized standard 
stars. The stars observed are listed in table [4] together with 
the parameters needed for the bandpass model. 

The wavelength dependence of interstellar polar- 
ization can be represented by t he empirical mod el 
llSerkowski. Mathewson fe Fordlll975l : IWilking et al.lfl980l ) 

P{ A) = Pma* exp (-K In 2 (19) 

\ Amax J 

where the values of Pmax, A max and K are empirical 
parameters that are fitted to observed wavelength dependent 
polarization measurements and are given in table [4] for our 
standard stars. 

The expected polarization in any of the HIPPI filters 
can then be predicted by averaging this function over the 
bandpass model as described in section El as follows: 


_ f P(\)e(\)S(\)d\ 
P f e(\)S(\)d\ 


( 20 ) 


Observations of polarized standard stars are listed in 
table [o] 

It can be seen from table [5] that the repeatability of 
these measurements from night to night is excellent with 
polarization agreeing to 0.01 % (100 ppm) or better and 
the position angles agreeing to typically 0.1 deg. The posi¬ 
tion angles also agree quite well with the expected values. 
While the position angle zero point was calibrated using 
these stars, the agreement between the three stars is an in¬ 
dicator of the accuracy of these measurements. The absolute 
calibration of position angle is somewhat poorer than the 
internal errors of ~0.1 degree, because it is limited by the 
available data on standard stars. 

Comparing with the predicted values of polarization 
from table[4]we find that the observed modulation efficiency 
averages ~90 % in the g’ band with somewhat lower values 
in the 500SP and r’ filters. These can be compared with the 
predicted efficiencies from the model described in section 
Eel which are given in the final column of the table. It can 
be seen that the observed efficiencies are in good agreement 
with the efficiencies predicted by our model. In most cases 
the measured and observed values agree within 1%, with the 
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Table 4. Predicted polarization in the HIPPI filters for polarized standard stars 


Star Spectral Type Eb—v Rv Pmax A max K 6 Expected P (%) Refs 

g' r' 500SP 


HD 23512 

A0 V 

0.34 

3.2 

2.29 

0.60 

1.06 

29.9 

2.15 


1, 2 

HD 147084 

A4 II/III 

0.72 

3.9 

4.34 

0.67 

1.15 

32.0 

3.77 4.26 

3.60 

3, 4 

HD 187929 

F6-G0 I 

0.18 

3.1 

1.76 

0.56 

1.15 

93.8 

1.70 

1.66 

5 


References: 1. iGuthrij <119871) , 2. iHsu &; Bregej dl982T) . 3. IWilking et alj <ll980l) , 4. iMartin, Clayton Sz Wol3 <1 19991) , 5. ISerkowski et all 

dl975h 


Table 5. Observations of polarized standard stars 


Star Date Filter Measured Expected Efficiency Efficiency 

P (ppm) 6 (deg) P (ppm) 8 (deg) (Observed) (Modelled) 


HD 23512 

Aug 28 

g' 

18852 

± 

37 

30.5 

± 

0.1 

21469 

29.9 

87.8 

90.0 

HD 147084 

Aug 29 

g' 

33919 

± 

11 

32.0 

± 

0.1 

37664 

32.0 

90.1 

91.0 


Aug 30 

g' 

33972 

± 

12 

32.0 

± 

0.1 

37665 

32.0 

90.2 

91.0 


Aug 30 

r' 

34910 

± 

17 

32.1 

± 

0.1 

42619 

32.0 

81.9 

84.2 


Aug 30 

500SP 

29349 

± 

17 

32.0 

± 

0.1 

35958 

32.0 

81.6 

80.7 


Aug 31 

500SP 

29337 

± 

12 

31.9 

± 

0.1 

35956 

32.0 

81.6 

80.7 

HD 187929 

Aug 28 

g' 

15450 

± 

9 

93.5 

± 

0.1 

16950 

93.8 

91.1 

91.2 


Aug 29 

g' 

15493 

± 

9 

93.6 

± 

0.1 

16960 

93.8 

91.3 

91.4 


Sep 2 

g' 

15560 

± 

13 

93.6 

± 

0.1 

16944 

93.8 

91.8 

91.1 


Sep 2 

500SP 

14213 

± 

16 

93.8 

± 

0.1 

16570 

93.8 

85.7 

84.7 


largest deviation being 2.3 %. Given that the quoted accu¬ 
racies of the standard star values are typically around 1% of 
the measured value this is excellent agreement. 


4.4 Polarization Sensitivity 

The polarization sensitivity of a polarimeter is its ability to 
measure small polarization levels, and is equivalent to the 
precision of measurement on low polarization objects. It can 
be estimated by looking at the night to night scatter of re¬ 
peat observations of low polarization stars. In table [6] we 
show such results for four sets of observations on three stars 
that have been measured on three or more nights. The data 
presented here are the instrumental Stokes parameters as 
measured by HIPPI before corrections for telescope polar¬ 
ization, and position angle zero point. 

The STDEV row at the bottom of the table gives the 
standard deviation of the numbers in the column and mea¬ 
sures that night-to-night scatter in the repeat observations. 
The average of these values is 4.0 ppm, equivalent to 4.3 
ppm when the modulation efficiency correction is included. 

While this represents our current best estimate of the 
sensitivity achieved with HIPPI, it should be noted that the 
internal errors of most of these observations are typically 
3-5 ppm, so part of the scatter may well be due to the 
contribution of random noise. For Beta Hydri (B = 3.4) the 
expected photon rate (above the atmosphere) at the AAT 
in the g' band is 1.06 x 10 10 photons sec -1 . Assuming a 
10 % total throughput allowing for atmospheric, telescope, 
and instrument losses as well as the detector QE (about 
25% averaged over the band), we detect a total of 1.92 x 
10 11 photons in the 180 second integration used for these 


observations. The expected photon shot noise limited sensi¬ 
tivity per Stokes parameter (1/ sjN p hotons) is then 2.3 ppm. 
The excess noise factor of the photomultiplier tubes which is 
typically around 1.2, will increase this to ~2.8 ppm. This is 
close to the measured errors for this object of 3.1 to 3.4 ppm 
indicating that even on these bright objects, where the pho¬ 
tomultiplier gain is set to its lower values, we are obtaining 
close to photon-noise-limited performance. 

The true polarization sensitivity achievable with HIPPI 
may therefore be somewhat better than these figures indi¬ 
cate. The lower figures for the brightest of these stars (Cano¬ 
pus) indicate a value nearer 3 ppm. 

These figures can be compared with the performance of 
other “parts- per-m illi on” p olarimeters. For PlanetPol, us¬ 
ing data from I Bailey et all d2008l ) the equivalent sensitivity 
fig ure is 2.1 ppm. For POLISH using data from table 3 of 
IWiktorowicz fe Matthewsl (120081) the figure is 8.5 ppm. The 
sensitivity of 4.3 ppm achieved with HIPPI shows that a po¬ 
larimeter using FLC modulation is quite competitive with 
the PEM technology used in these other two instruments. 


5 CONCLUSIONS 

We have built and tested a stellar polarimeter based on 
a ferro-electric liquid crystal (FLC) modulator. The po¬ 
larimeter has been used successfully on the Anglo-Australian 
Telescope. The FLC does not provide as “polarimetrically 
clean” a system as the photo-elastic modulators (PEMs) 
used in other high precision polarimeters. Nevertheless, af¬ 
ter correction of instrumental effects using multiple stages of 
modulation, HIPPI achieves performan ce comparable with 
PEM polarimeters such as PlanetPol (|Hough et al.l 120061 ) 
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Table 6. Repeat observations of low polarization objects to determine sensitivity 


Star Canopus (g') Beta Hyi (g') BS 5854 (g') Beta Hyi (500SP) 

Q/I (PPm) U/I (ppm) Q/I (ppm) U/I (ppm) Q/I (ppm) U/I (ppm) Q/I (ppm) U/I (ppm) 


130.7 ± 0.8 

-62.7 ± 0.8 

42.1 ± 3.2 

-22.7 ± 3.1 

28.5 ± 4.8 

-14.5 ± 2.4 

38.5 ± 3.4 

-27.2 ± 3.4 

134.4 ± 1.2 

-64.9 ± 1.8 

42.2 ± 3.2 

-37.8 ± 3.2 

37.7 ± 3.6 

-17.5 ± 13.8 

41.0 ± 3.3 

-27.6 ± 3.3 

129.7 ± 1.3 

-58.4 ± 1.3 

43.1 ± 3.3 

-31.3 ± 3.1 

41.2 ± 3.6 

-19.9 ± 3.3 

35.2 ± 3.5 

-19.6 ± 3.6 

128.7 ± 3.1 

-58.7 ± 2.5 

35.0 ± 3.4 

-32.4 ± 3.1 

33.5 ± 5.1 

-23.3 ± 3.6 



STDEV 2.5 

3.2 

3.8 

6.2 

5.5 

3.7 

2.9 

4.5 


and POLISH (Iwiktorowicz fc Matthews! 12008 ). while offer¬ 
ing advantages of lower cost, compact size and greater effi¬ 
ciency. 

HIPPI provides both high precision - it can mea¬ 
sure the polarization of highly polarized stars to a precision 
in fractional polarization of 0.01 % and position angle to 
around 0.1 degree or better — and high sensitivity — it can 
measure very low levels of polarization with a precision on 
low polarization objects of 4.3 ppm or better. The wide dy¬ 
namic range of its detectors allow it to observe objects from 
the brightest stars in the sky (e.g. the observations of Sirius 
and Canopus reported here), while still having the capabil¬ 
ity to observe quite faint objects (dark noise from the PMTs 
should only start to limit performance for stars fainter than 
about 16th magnitude). 

We have measured the telescope polarization of the 
AAT to be 48 ± 5 ppm in the SDSS g' filter. The preci¬ 
sion of HIPPI is such that the primary limitation on the 
absolute accuracy of its measurements is currently set by 
the lack of suitable standard stars for calibration. There are 
very few low polarization stars measured at the parts-per- 
million level, and this limits our ability to reliably determine 
and correct for the telescope polarization. Similarly the pre¬ 
cision of our measurements of polarized standard stars is 
better than the absolute calibration of their polarizations 
and position angles currently available. 
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